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Gaseous X-ray polarimetry refers to a class of detectors used for measuring the polarization of soft X-rays. 
The systematic effects of such detectors introduce residual modulation, leading to systematic biases in the 
polarization detection results of the source. This paper discusses the systematic effects and their calibration and 
correction using the Gas Microchannel Plate—Pixel Detector(GMPD) prototype for Polar2/Low Energy X-ray 
Polarization Detector(LPD). Additionally, we proposes an algorithm that combines parameterization with Monte 
Carlo simulation and Bayesian iteration to eliminate residual modulation. The residual modulation after data 
correction at different energy points has been reduced to below 1%, and a good linear relationship is observed 
between the polarization degree and modulation degree. The improvement in modulation degree after correction 
ranges from 2% to 15%, and the results exceed those of IXPE above 5 keV. 
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I. INTRODUCTION 


In recent decades, there has been a growing interest in 
the field of gamma-ray astronomy, particularly in the study 
of Gamma-ray bursts (GRBs)[1], which are highly energetic 
cosmic events. While satellite observations from missions 
like Swift[2] and Fermi[3] have yielded valuable insights into 
the energy spectra and timing properties of GRBs, numerous 
fundamental questions remain unanswered. These include 
understanding the mechanisms that propel the energetic jets, 
elucidating the processes responsible for energy dissipation, 
determining the composition of the jets, investigating the con- 
figurations of magnetic fields, and unraveling the mechanisms 
behind particle acceleration and radiation[4—9]. The detec- 
tion of polarization in GRBs plays a crucial role in providing 
important clues for addressing the aforementioned issues[10— 
14]. 

Scheduled for deployment in 2026 as an external payload 
on the China Space Station, POLAR-2[15] is the successor 
experiment to POLAR[16]. Its main goal is to conduct high- 
precision measurements of polarization across the spectrum 
from soft X-rays to gamma rays. The GMPD[17] is an inno- 
vative gas pixel detector developed to validate the design of 
the POLAR2/LPD[18] payload. 

Recently launched polarimetric detectors such as 
PolarLight[19], IXPE[20], and the under-development 
eXTP[21] and Polar-2/LPD all utilize gas pixel polarimetric 
detector structures. This type of detector has high spatial 
resolution, capable of imaging electron tracks at the level of 
hundreds of micrometers, thus providing excellent sensitivity 
in polarimetric detection. However, due to the complex and 
sophisticated structure of the gas pixel detector, as well as its 
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high spatial resolution sensitivity, the operational state of the 
instrument, various components, and electronic devices may 
introduce systematic effects on the polarimetric detection 
results. These systematic effects can result in non-zero 
modulation named residual modulation when detecting 
unpolarized sources, leading to systematic biases in the mea- 
surement of polarized sources. Since low-energy electron 
tracks are relatively short, the residual modulation effects 
produced by these systematic effects are more significant in 
low-energy events and cannot be ignored. 


For residual modulation, IXPE employs two methods for 
correction[22]: the first involves oscillating the detector dur- 
ing the detection process to integrate and eliminate some of 
the systematic effects. The second method involves calibrat- 
ing the corresponding Stokes parameters q and u for sys- 
tematic effects in different regions and energy points, and 
then subtracting them on a event-by-event basis to eliminate 
the systematic effects. In this reasearch, We listed some of 
the known causes of residual modulation and corrected some 
of these effects based on their generation mechanisms. For 
another part of the residual modulation where the specific 
causes are currently unclear, we proposed a correction algo- 
rithm and obtained favorable outcomes. 


In this paper, we first introduced the basic structure and 
polarization detection principles of the Polar2/LPD detector 
in Sect.II. We then discussed the residual modulation caused 
by signal response and its correction methods in Sect.II. 
In Sect.IV, we discussed the residual modulation caused 
by geometric effects and proposed a modulation curve cor- 
rection method based on the parameterization of scale ra- 
tios, combined with Monte Carlo simulation and Bayesian 
iteration[23](see in Appendix A), and provided the errors of 
this algorithm. Subsequently, we compared various data re- 
construction characteristics before and after algorithm cor- 
rection and compared them with the modulation calibrated 
by the IXPE detector. Finally, in Sect.V, we discussed the 
performance of the GMPD after correction, emphasizing the 
performance and scalability of the correction algorithm, and 
outlined prospects for future work. 
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II. GEOMETRIC STRUCTURE AND WORKING 


PRINCIPLE OF LPD 


The LPD system shown in Fig.l is composed of a total 
of 9 detector modules, arranged in a 3x3 array configura- 
tion. Each detector module consists of 9 detection units with 
90° field of view (FoV), resulting in a total of 81 detection 
units. A detection unit of LPD consists of working gas, a Gas 
Micro-Channel Plate (GMCP)[24], a pixel readout chip, and 
frame structure. Working gas serving the purpose of photo- 
electric effects and formation of ionization tracks. The up- 
per end of the gas chamber is sealed with a 50 um beryllium 
window, which prevents entry of lower energy photons and 
ensures gas containment to prevent leakage. GMCP layer is 
positioned near the bottom plane at the chamber, for the pur- 
pose of electron avalanche multiplication. At the bottom of 
the chamber is the chip. The intended chip version, Topmetal- 
L, which is expected to be formally integrated into LPD, will 
undergo optimization in terms of power consumption, effec- 
tive area, and resolution based on the existing chip version, 
Topmetal-H[25]. 

The soft X-rays, as shown in Fig.2, pass through the beryl- 
lium window of the detector unit and have a certain probabil- 
ity of undergoing photoelectric effects within the drift region, 
resulting in the generation of photoelectrons. These photo- 
electrons carry the polarization information of the incident 
photons. Photoelectrons deposit ionization energy within the 
gas and generate secondary ionization electrons until they 
come to a complete stop. Within the induction region, an 
upward-directed electric field is applied, causeing some of 
the secondary ionization electrons drift downward onto the 
surface of the GMCP. A portion of these electrons enters the 
micro-channels and undergoes avalanche multiplication. The 
multiplied electrons then emerge from the lower surface of 
the GMCP, where some of them are absorbed, resulting in 
the production of a pulse signal. The remaining multiplied 
electrons continue to drift towards the Topmetal chip, induc- 
ing signals in the corresponding pixel positions. This process 
allows the projection of the photoelectron track onto the 2D 
plane of the Topmetal chip, enabling us to obtain the projec- 
tion results. 

In general, the angular distribution of photoelectrons de- 
tected by the gas pixel detector is modulated by polarized 
X-rays. For gas pixel detectors, photons primarily interact 
with the K-shell electrons of gas molecules through photo- 
electric interactions, and the direction of electron emission 
is described by the differential cross-section according to the 
following formula[26]: 
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where ( is the emission velocity of the photoelectron in 
units of the speed of light c, 0 and ¢ are the latitude and az- 
imuth angles, respectively. Due to the lack of resolution in the 
Z direction for the LPD detector, we reconstruct the 2D pro- 
jection angular distribution of the azimuthal angle of the pho- 
toelectron emission, corresponding to the integration of 0 in 
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the formula. Therefore, the reconstructed angular distribution 
is modulated by the cos? factor. In theory, for 100% polarized 
X-rays, the minimum value of the true emission distribution 
in phi should be 0. However, due to limitations in instru- 
ment resolution, system effects, and reconstruction algorithm 
accuracy, there will be a certain proportion of unmodulated 
components in the angular distribution, as shown in Fig.3 (a), 
(b). Therefore, the modulation function of ¢, M(¢) can be 
written as: 

M(¢) = A + Beos(2(¢ — ¢o)). (2) 
Therefore modulation factor u is defined as the ratio of the 
area occupied by the modulation component in the distribu- 
tion: 
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area, (3) 


u 


HI. CALIBRATION AND CORRECTION OF SIGNAL 


RESPONSE 


The structural design and operational principles of GMPD 
result in variations in the response between pixels, which can 
impact the energy resolution of the detector. More impor- 
tantly, some of these factors can introduce anisotropic differ- 
ences, leading to residual modulation. This section primarily 
discusses the impact and calibration of these factors. 


A. Pixel response differences 


Due to the subtle structural differences between each pixel, 
the uniformity of the electric field, and the uniformity of 
GMCP gain, the signal induction intensity of drift charge 
varies among different pixels. It is necessary o calibrate the 
relatives ignal induction intensity on each pixel. We uni- 
formly irradiate with a 4.51 keV flat source and statistically 
record the signal distribution received by each pixel. As the 
response curve of the pixels exhibits good linearity[27], we 
can characterize the relative signal induction intensity of a 
pixel by the mean of the signal distribution received on that 
pixel. In the process of calculating the mean, we only selected 
the part of the signal intensity greater than 50 in order to elim- 
inate the interference of noise. Fig.4 illustrates the average 
distribution of pixel ADC values before and after correction. 


B. Rolling-Shutter and Signal Decay 


Another source of residual modulation is the attenuation of 
pixel signal amplitude caused by signal readout time delay. 
Since Topmetal-II adopts a rolling-shutter readout method to 
read each frame of the image, pixel signals are read out in 
sequence, which means there is a certain delay from the trig- 
gering of the induction signal to the readout, and there is also 
a delay in the readout time of different pixels on the same 
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Fig. 4. The distribution of pixel average ADC values before and after 
correction, with blue representing the pre-correction values and red 
representing the post-correction values. 
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of scanning the chip column by column along the 0° direc- 
tion, the scanning time interval between adjacent pixels in 
the 90° direction is 35 us, while the scanning time interval 
between adjacent pixels in the 0° direction is 500 ns, with 
AToo < AToo° . The difference in scanning time intervals 
between 0° and 90° can result in inconsistent signal attenua- 
tion gradients in these two directions, introducing a vertical 
bias, namely, residual modulation in the 90° direction. 

In order to calibrate the systematic errors caused during the 
scanning process, it is first necessary to calibrate the signal 
attenuation behavior of each pixel, and secondly to determine 
the time difference between each triggered and readout pixel. 
We input square wave signals to the chip and record the output 
results of the pixel readout signals for multiple consecutive 
frames, in order to obtain the decay characteristics of each 
pixel and perform parameter fitting. The theoretical formula 
for pixel decay is given by Equation 4: 


ae 
b-t+ec) 


Fig.5(b) illustrates the decay pattern of signal intensity over 
time on a pixel and the fitting result. Fig.5(c) shows the decay 
time distribution of all pixels on Topmetal-II, indicating that 
the typical decay time scale for pixels is 20 ms. 

The time precision of GMCP is 10 ns [27]. The time differ- 
ence At is obtained by comparing the trigger signal of GMCP 
and the trigger signal position on Topmetal-II. Since the typ- 
ical time scale required for electron propagation in GMPD is 
on the order of tens of ns, it can be neglected compared to 
the characteristic time scale of pixel decay. Based on the time 
difference, we can then correct the decay signal for each pixel 
using the following formulas: 


A(t) =a- exp (- (4) 


to — At 
Atruth = @- : 
ai exp ( TORET] a 
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0 1—b-log(42-) (6) 


C. Charge pile up effect 


The encapsulated detector exhibits an initial stage where 
the gain increases with the accumulated number of events, as 
illustrated in Fig.6(a). This effect is attributed to charge ac- 
cumulation. The surface of the Topmetal-II utilized in the de- 
tector features a grid-like insulating layer, causing electrons 
that fall and become adsorbed on this layer to have limited 
mobility. As the accumulation of avalanche multiplied elec- 
trons rises, the potential on the chip surface gradually alters, 
impacting the charge collection efficiency and modifying the 
gain, as depicted in Fig.6(b). 

If the charge accumulation process is unevenly distributed 
on the surface of the chip, it will lead to a noticeable sig- 
nal intensity gradient on the chip surface, eventually result- 
ing in the generation of pseudo-modulation perpendicular to 
the gradient direction.Fig.7 illustrates the residual modulation 
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caused by the charge accumulation effect. Initially, a ferrous 
strip was used to partially obstruct a section of the detector’s 
field of view, leaving a gap of a few millimeters. Follow- 
ing a 2-hour exposure to an X-ray flat source, the obstruc- 
tion was removed, and a 5.9keV unpolarized Feř source was 
used to irradiate and collect the photoelectron tracks. Upon 
reconstruction, it was observed in Fig.7(a) that the signal gain 
at the previous narrow gap position was significantly higher 
than the shaded area, and the residual modulation in the nar- 
row gap area was higher than in the shaded area, with the 
modulation direction parallel to the gap. Subsequently, with- 
out any obstruction, the X-ray flat source was used again for 4 
hours to accumulate charges on the entire surface of the chip 
to near saturation. The detector was then irradiated with the 
5.9keV unpolarized Fe®®, and the tracks were reconstructed in 
Fig.7(b). Comparing the results of the Fe55 measurements be- 
fore and after charge accumulation reached saturation, it was 
found that the residual modulation caused by the uneven gain 
due to charge accumulation significantly decreased. There- 
fore, it is possible to mitigate the impact of the charge accu- 
mulation effect by calibrating or measuring the detector after 
saturating the charge accumulation before conducting exper- 
iments. Since the accumulated charge is unlikely to naturally 
dissipate, once the detector is encapsulated, only one thor- 
ough charge accumulation is required. 


By employing Garfield++ and COMSOL for charge drift 
accumulation iteration and updating of the drift electric field, 
we successfully replicated this effect in simulations, as indi- 
cated by the blue data points in Fig.6(a), which align with the 
experimentally observed gain variation results. The process 
of charge accumulation can be described by a simplified Eq.7. 
Where n is the number of events, q is the accumulated charge 
on the chip, and qmax is the maximum saturated accumulated 
charge, and a, is the charge adsorption coefficient. Therefore 
the change in the accumulated charge quantity with respect 
to the detector counts, q(n), can be expressed in a parametric 
form as given in Eq.8. 


_ =o (1 z a) l (7) 
n max 
q(n) = zo + z1 exp(#o(n + z3)). (8) 


Both experimental and simulation results indicate that the 
charge accumulation process gradually reaches saturation, 
leading to a stable final gain. Additionally, the non-focusing 
observation mode of the LPD can prevent the uneven accu- 
mulation of charge on the chip surface. Therefore, after a 
sufficient number of accumulated events, the impact of the 
pseudo-modulation caused by the charge accumulation effect 
in the encapsulated detection unit can be reduced to a negli- 
gible level. 
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Fig. 5. (a) Topmetal-II~ pixel CSA structure. The green portion represents the top metal, which is used to induce the drift charge signal. Cr 
denotes the feedback capacitor. V, represents the gate voltage, Va represents the drain voltage, and Vef represents the amplifier’s reference 
voltage. (b) Trend and fitting curve of pixel ADC value decay over time (c) Distribution of decay times for all pixels, where the decay time of 
a pixel is defined as the time interval for the pixel signal value to decay from a to a/2. 


14000 


13000 4 


12000 4 


11000 4 


10000 4 


Peak ADC 


9000 4 


8000 + 


— Fit 
COMSOL & Garfield++ Simu 
@ Topmetal Peak ADC 


7000 4 


6000 + 


60000 80000 100000 120000 


Counts 


(a) 


0 20000 40000 


Topmetal-ll 
(b) 


Fig. 6. (a) The variation of chip collection efficiency with the accumulation of events. The orange points represent experimental data results 
of 4.51 keV photoelectron track energy deposition, characterized by fitting the peak ADC of the data spectrum to represent the variation in 
chip charge collection efficiency. The blue points represent the trend of collection efficiency variation obtained from the joint simulation using 
COMSOL and GARFIELD++. (b) Charge Accumulation Schematic: Electrons accumulate on the insulating layer, forming a low potential 
region around the top metal. This leads to the formation of a funnel-shaped electric field above the top metal, thereby enhancing the charge 


collection efficiency of the pixel. 


IV. CALIBRATION AND CORRECTION OF 
GEOMETRICAL EFFECTS 


A. Pixelization influence 


As shown in Fig.8(a), we consider a shorter track with 
a circular projection. Due to the parallel arrangement of 
Topmetal-II chips in the X and Y directions, the signal distri- 
bution sensed on the chip pixels exhibits anisotropy for such 
tracks. The symmetry is most pronounced in the directions of 
0° and 90°, which are aligned with the pixel arrangement. The 
commonly used moment analysis algorithm for such shorter 
tracks calculates the centroid line of the pixel track to deter- 
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mine the direction of electron emission. This can lead to a 
bias in the reconstruction direction of these tracks towards 0° 
and 90°. 


To mitigate the residual modulation caused by pixel ar- 
rangement, we need to exclude events with too few respon- 
sive pixels and events that are too short or too circular during 
event selection. Therefore, during reconstruction, we only se- 
lect events with a number of hit pixels greater than or equal 
to 27 and exclude the bottom 20% of events with smaller el- 
lipticities. Fig.8(b) below shows the angular distribution of 
the reconstructed unbiased events before and after the event 
selection. After the event selection, the residual modulation 
caused by pixel arrangement is significantly improved. 
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Fig. 7. (a) The residual modulation distribution of 5.9keV Fe°° tracks after uneven charge accumulation due to narrow gap obstruction. 
(b) The residual modulation distribution of 5.9keV Fe°° tracks after uniform charge accumulation following the removal of the narrow gap 
obstruction. The heatmap represents the distribution of the reconstructed photoelectron emission positions of the tracks, with the direction of 
the red lines indicating the direction of residual modulation. The length of the line and the adjacent number represent the value of the residual 


modulation. 
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Fig. 8. (a)The shorter circular tracks are projected onto the read- 
out Topmetal-II chip, where the pixel grayscale represents the signal 
intensity, with darker colors indicating stronger signals. (b)The an- 
gular distribution of reconstructed electron tracks at 2.98keV, with 
a gas environment consisting of a 40% volume fraction of helium 
gas and 60% dimethyl ether (DME) at 0.7atm. The red data points 
represent the reconstruction results without any filtering based on 
the number of triggered pixels and ellipticity. The green data points 
represent the filtered results, excluding events with a number of trig- 
gered pixels less than 36 and the bottom 20% of events with smaller 
ellipticity. 


B. Truncation effect 


Similarly, due to the Rolling-Shutter line-by-line scanning 
readout of the chip, if an event occurs precisely at the po- 
sition covered by the pixels being scanned at that moment, 
the event will be truncated and appear in both the preceding 
and subsequent frames. If the truncated part in one frame has 
fewer fired pixels that do not exceed the threshold, we can 
only obtain an incomplete truncated event. As the edge of the 
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truncated track is always parallel to the scanning direction, it 
introduces a systematic bias in scan direction. 

Thus, we continue to utilize the time information from 
GMCP and Topmetal-II to determine if an instance is trun- 
cated. Considering the combined time resolution of GMCP 
and Topmetal-II is 262 ns[27] and Tpizeļ, we determine if a 
track is truncated by examining whether the pixel scanned 
when the signal arrives and the positions of the five pixels be- 
fore and after it precisely overlap with the region covered by 
the photoelectron track signal. 


C. Track image distortion 


Excluding the systematic effects and corrections discussed 
above, the angular reconstruction of track data obtained from 
the detector still exhibits some residual modulation. This may 
partly be attributed to the geometry and potential distribution 
of the detector. The gas cavity of the LPD detection unit is not 
completely symmetrical. As shown in Fig.9, in addition to the 
charge induction chip Topmetal, a temperature and pressure 
sensor chip is also placed nearby. This placement leads to 
a relatively significant distortion of the electric field near the 
side of the Topmetal chip adjacent to the sensor chip, result- 
ing in a noticeably higher residual modulation on that side. 
Furthermore, there is a 1 mm wide and 0.8 mm deep groove 
between the charge induction collection plane of the Topmetal 
chip and the base plane of the detection unit. Additionally, 
several to a dozen bonding wires are present around the chip. 
The geometric structure of the chip’s edge and the potential 
on the bonding wires also cause distortion of the electric field 
at the edge of the chip. Consequently, it can be observed that 
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the direction of the residual modulation reconstructed in the 
edge portion of Fig.7 is generally perpendicular to the edge 
of the chip. Therefore, in order to minimize the influence of 
edge electric field distortion on the reconstruction, we choose 
to exclude events within 12 pixels of the charge center dis- 
tance from the edge when selecting valid events. 


Fig. 9. (a) Physical diagram of the detector unit base, with the larger 
central chip being the Topmetal-II charge induction chip, and the 
smaller chip on the left being the temperature and pressure sensing 
chip. (b) Simulation results of the electric field near the detector 
base. Due to the structural characteristics of the pressure sensing 
chip, the overall geometric structure of the detection unit exhibits a 
certain degree of asymmetry, leading to some distortion in the elec- 
tric field. 


The residual modulation distribution in different regions 
near the center of the chip appears to be more random. This 
variability in residual modulation in certain regions may stem 
from systematic process issues during chip etching, subtle ir- 
regularities during detector installation, and the uneven accu- 
mulation of charge resulting in differences in the electric field 
across different areas of the chip. These issues can all impact 
the electric field distribution near the chip surface, and the 
distortion of the electric field can alter the track shape. This 
alteration is often nonlinear, and the impact on tracks at dif- 
ferent positions, heights, and lengths varies. As a result, we 
lack sufficiently precise information to make pixel-by-pixel 
or event-by-event corrections from a first-principles perspec- 
tive for the obtained tracks in the experiment. 

The deformation of tracks is reflected in the differences in 
position resolution in different directions of the detector. As 
shown in the Fig.10, at different energy points, the position 
resolution in the X direction of the detector is worse than that 
in the Y direction. This indicates that the distortion of the 
track is more severe in the X direction, and this anisotropic 
deformation of tracks leads to excessive stretching in the X 
direction, resulting in significant residual modulation. 

Similar residual modulation issues also arise in the IXPE 
detector. The correction scheme for residual modulation pro- 
vided by IXPE involves calibrating experimental scales for 
each chip region to correct the Stokes parameters required for 
event reconstruction. Since the IXPE detector needs to im- 
age the observed objects, segmenting and correcting differ- 
ent regions is necessary. However, for the LPD, which lacks 
imaging capabilities, photons from the source will uniformly 
fall on the entire chip surface. Therefore, the LPD only needs 
to consider correcting the distribution of residual modulation 
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Fig. 10. The position resolution of GMPD of varying energies. Cir- 
cular markers depict the results in the X-direction, while triangular 
markers represent the results in the Y-direction. 


integrated over the entire chip surface for events. To address 
this, we propose a Bayesian method combined with Monte 
Carlo simulations to correct residual modulation. 


1. Correction algorithm 


When correcting the data for an energy point, we only need 
to calibrate a correction parameter 7): the ratio of the pixel size 
in the Y direction to the pixel size in the X direction. We can 
phenomenologically explain the need to introduce the param- 
eter 77: the distortion of the electric field will cause the equipo- 
tential surfaces to no longer be parallel to the Topmetal chip 
plane. Therefore, by projecting the chip plane onto the de- 
formed equipotential surface, the linearity in different direc- 
tions of the chip will have different scaling rates. We select 
the ratio of the scaling rates calibrated in the two directions 
parallel and perpendicular to the scanning direction as 7. It 
should be noted that the 7) value for different regions of the 
chip is different. However, because the LPD does not have 
polarized imaging capabilities, the correction parameter we 
consider is actually the weighted average value 77 of the pa- 
rameters for different chip regions. 

Calibrating 7 requires a comparison of experimental data 
and simulations. Taking the 5.40 keV energy point as an ex- 
ample, we calibrated the modulation curves obtained from a 
99.9% polarized source at 0°, 30°, 60°, 90°, 120°, and 150° 
polarization phases. It can be observed that due to residual 
modulation, there are significant differences in the modula- 
tion at these phases, with a difference of approximately 18% 
between the modulation at 0° and 90° as shown in Fig. 15(g). 
Since residual modulation is an inherent property of the de- 
tector and is independent of the polarization phase of the 
source, the overall modulation curve measured is a result of 
the superposition of residual modulation and source modu- 
lation. Therefore, the overall modulation curve can be de- 
scribed as: 
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Movs(¢) = Mres(¢) . Msource( Q, ġo). (9) 


Where Mops is the modulation curve obtained from recon- 
structed angular distribution data, res represents the impact 
of residual modulation, and Msource is the modulation curve 
generated by the polarized source. At normal incidence, the 
form of Mspource 18: 


Mgource(¢, ġo) = Acos (b E po)” +B. 


According to equation 10, we observe that Msource is mod- 
ulated by cos?. Therefore, by equally combining two sets of 
data with a 90° difference in polarization phase, the modu- 
lation caused by the polarized source can be eliminated. As 
a result, when equally mixing six sets of data at 0°, 30°, 60°, 
90°, 120°, and 150° polarization phases, the modulation curve 
of the angular distribution Mops « Mpres. The modulation 
distribution of the combined data is shown in Fig.11. The 
combined results indicate that the residual modulation dis- 
tribution still follows equation 10, and fitting different com- 
bined data sets within the error range shows that the residual 
modulation values obtained from different data sets are con- 
sistent, with the phase of the residual modulation being 0°. 
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Fig. 11. Combination results of data at different polarization phases 
for 5.40 keV. 


We obtained the residual modulation amplitude at the 
5.40 keV energy point by fitting the residual modulation curve 
from Fig.11(d). Next, we consider using simulations to re- 
produce the same residual modulation distribution and ob- 
tain a response matrix for correcting the residual modulation 
in the experimental data. We utilized the star-XP software 
framework[28] specifically designed for the LPD detector. 
The simulation framework meticulously simulates the inter- 
action processes between photoelectrons and the detector, as 
well as the digitization process. The simulated data output by 
the framework shows good agreement with the experimental 
data. Our operational procedure followed the steps below: 
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1. In the simulation framework, we simulated the tracks 
of 1.5 million unpolarized 5.40 keV X-ray photons and 
maintained the parameters set in the simulator consis- 
tent with the actual operating parameters of the detec- 
tor. 


2. Initially, we set 7=1, representing the state of the detec- 
tor without electric field distortion, and simulated the 
two-dimensional image of the photoelectron tracks af- 
ter digitization. We reconstructed each of the 1,500,000 
tracks without distortion to obtain the reconstructed an- 
gle information, Anglen. It is important to note that 
this is not the true value of the emission angle of the 
photoelectrons provided by the simulation, but rather 
the angle value obtained from the reconstruction. The 
distribution of Angler is shown in Fig. 12(a). 


3. In the simulation, we adjusted the value of 7) to deviate 
from 1, representing the occurrence of electric field dis- 
tortion in the detector. We used the photoelectron track 
simulations from Step 1 after digitization, and due to 
different scaling in the X and Y directions, the recon- 
structed angle distribution, Anglep,,,,,, exhibited a non- 
zero residual modulation. When 7j<1, the phase of the 
residual modulation is 0°, consistent with the experi- 
mental data. By adjusting the value of 7, we were able 
to align the modulation amplitude of the Anglepistor dis- 
tribution with Fig.11(d), as shown in Fig.12(b). For 
5.40 keV, the value of 77 was determined to be 0.981. 


. Combining the Angle,,,,,, and Anglepistor reconstructed 
step by step in the second and third steps, we can ob- 
tain the response matrix Mpistor, which arises due to the 
adjustment of the parameter 7. The physical interpreta- 
tion of Mpistor is as follows: if we denote i as the 
number of instances for which the reconstructed angle 
falls within the i-th bin when 7=1, then the number of 
instances for the same events after adjusting 7 and re- 
constructed within the j-th bin is given by equation 11. 
In other words, the element (m,n) of Mpistor is propor- 
tional to the probability value P(nlm): the probability 
that an event originally reconstructed in the m-th bin is 
reconstructed in the n-th bin due to the adjustment of 77. 


(1) 


After obtaining the parameter 7 1 that describes the 
system effect and its corresponding response matrix Mpistor, 
we employed a Bayesian iterative process algorithm to de- 
couple the modulation distribution generated by the polar- 
ization source, which does not have electric field distor- 
tion, from the overall system effects. Many software pack- 
ages offer Bayesian algorithm capabilities, and in our study, 
we utilized the RooUnfoldBayes packages integrated within 
RooUnfold[29]. The RooUnfoldBayes package is capable of 
iteratively restoring the input angular modulation distribution 
to its undistorted state, based on the provided response ma- 
trix Mpistor, and automatically calculates the errors for each 
bin following the Bayesian iteration process. 
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Fig. 12. (a) The distribution of Angley,,,,, with the red line representing the fitted modulation curve, exhibiting a modulation degree of 0. (b) 
The distribution of Angle,,.,,,. With the modulation degree adjusted by tuning 77 to match the experimental data in Fig.11(d). (c) The response 
matrix, with the ordinate representing the reconstructed angles of events at 7=1, and the abscissa representing the reconstructed angles of 


events at 7=0.981. 


The use of Bayesian methods involves the selection of prior 
distributions and the adjustment of the number of iterations. 
Firstly, due to the periodicity of the modulation curves with a 
4 period of 7, monotonically increasing or decreasing distribu- 
tions are not appropriate. Therefore, for simplicity, we set the 
prior distributions to uniform distributions. Secondly, con- 
cerning the number of iterations, we determine the conver- 
gence of the iteration process by comparing the x? values of 
the distributions M(@)n+1 and M(¢@), after the n + 1 th and 
n th iterations. We found that when the number of iterations is 
set to 10, the x? values for different phases, polarizations, and 
energies are all less than 0.5, which indicates that the iterative 
process has essentially reached convergence. Additionally, 
after 10 iterations, the introduced iteration errors in each bin 
are relatively small. Therefore, we set the number of itera- 
tions to 10. Fig. 13(a) illustrates the variation of the x? values 
corresponding to different numbers of iterations, while Fig.s 
13(b) and 13(c) present the corrected results for the 5.40 keV 
99.9% polarized data at 0° and 90° phases, respectively. 
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sis At a specific energy point, using the aforementioned 
sis method, we only need to calibrate one corresponding param- 
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Ey 


eter, namely 7, and simulate the response matrix at that en- 
ergy point. This allows for the application of Bayesian itera- 
tion to correct the modulation curves at different polarization 
degrees and phases. The corrected polarization degree and 
modulation level exhibit a good linear relationship, and the 
modulation levels at different phases also show good consis- 
tency. Fig.14 illustrates the comparison of the modulation 
distribution before and after correction for experimental data 
at 5.40 keV and 90° polarization phase, ranging from unpolar- 
ized to 99.9% polarized. When the polarization degree is low, 
the residual modulation will dominate the distribution of the 
modulation curve. At this point, the unrevised experimen- 
tal data reconstruction results will exhibit significant devia- 
tions. In contrast, the corrected data maintains good stability 
in the reconstruction of the polarization phase, while also ex- 
hibiting a strong linear relationship between the modulation 
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s533 level and the polarization degree. Fig.15 displays the com- 
s534 parison of experimental data before and after correction at 
several energy points, including 2.98 keV, 4.51 keV, 5.40 keV, 
6.40 keV, and 8.05 keV, for different polarization phases and 
degrees. The uncorrected data, due to the residual modula- 
tion not being eliminated, exhibit significant differences in 
the reconstructed modulation degree at different phases. For 
example, at 4.51 keV, the modulation degrees at 0° and 90° 
for the same fully polarized source differ by approximately 
18%. The polarization phase reconstruction results from the 
uncorrected data also show a significant deviation from the 
true values, especially in the direction that differs by 90° from 
the polarization direction of the residual modulation. How- 
ever, after Bayesian iterative correction, the polarization data 
show good consistency in modulation degree across different 
phases. Additionally, the polarization degree and modulation 
degree exhibit a good proportional relationship, meeting the 
calibration requirements of the LPD. 
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The error in the modulation degree of the corrected data 
distribution mainly arises from two sources. One part origi- 
nates from the statistical error of the data, which can be ob- 
tained through fitting. The other part of the error arises from 
the process of using the Bayesian method for correction: 
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1. Error propagation in the Bayesian iteration process: 
This error can be calculated through the error propa- 
gation matrix A4. 
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. Termination of the Bayesian iteration: Although the 
chi-square calculation results show good convergence 
after 10 iterations for all experimental data, the conver- 
gence levels of the data at different polarization phases 
are inconsistent due to the fixed number of iterations. 
This results in slight differences in the reconstructed 
modulation degree at different polarization phases af- 
ter correction. 
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Fig. 13. (a) The x? variation of the distributions after 1-10 Bayesian iterations for different polarization degrees and polarization phases 
at different energy points. After 10 iterations, the x? values are all less than 0.5, indicating that the Bayesian iterations have essentially 
converged. (b) and (c) respectively show the completely polarized data at 5.40 keV and 0°, 90° phases. The green curve represents the original 
measured reconstructed angular distribution, while the blue curve represents the distribution after 10 Bayesian iterations. A comparison 
between (b) and (c) reveals that the Bayesian iterations have corrected the modulation levels at the two phases to the same level. 
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Fig. 14. (a)-(k) Comparison of modulation degree before and after correction for experimental data at 5.40 keV and 90° polarization phase 
for different polarization degrees where PD represents the polarization degree. The unpolarized data were obtained by mixing datasets with a 
90° phase difference using the method described at Fig.11, and the datasets for different polarization degrees were obtained by proportionally 
mixing unpolarized and fully polarized data. (1) The triangle represents the relationship between the modulation degree reconstructed from 
experimental data before correction and the polarization degree, while the circle represents the relationship after correction, and the red line 
represents the fitted curve of the corrected data points. 


3. Parameterized response matrix: The error in the esti- s4 The error propagation in point 1 is calculated by the 
mation of the parameter 7 provided by the simulator sss RooUnfold package. For the statistical error of the data and 
will be propagated to the response matrix, and during 576 points 1 and 2, due to the dependence of the Bayesian method 
the Bayesian iteration process using the response ma- s77 iteration process on the original data itself, it is difficult to 
trix, the error will be propagated to the corrected data’s s7s decouple and analyze the contributions of these two parts. 
modulation degree. s7 Therefore, a unified error estimation is provided using the 
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Fig. 15. The results are for 2.98 keV, 4.51 keV, 5.40 keV, 6.40 keV, and 8.05 keV, respectively. (a), (d), (g), (j), (m) represent the modulation 
degree measured at different polarization phases of polarized sources and the corrected modulation degree. (b), (e), (h), (k), (n) represent the 
comparison before and after correction for data at different polarization phases and degrees. The hollow triangles represent the uncorrected 
results, the hollow circles represent the corrected results, with different colors representing different polarization phases, and the red line 
represents the linear fitting of the corrected data. (c), (f), (i), (1), (0) represent the reconstruction of the polarization phase before and after 
correction, with the dashed line representing the true polarization phase of the polarized source. 
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modified Bootstrap method, and this part of the error is de- 
noted as Ounfola: Sampling 10,000 times at a certain polariza- 
tion degree (taking fully polarized data as an example). 


1. Each sampling involves 100,000 with-replacement 
samplings of the data at 0°, 30°, 60°, 90°, 120°, and 
150° phases in the experiment. 


2. Reconstruction of the sampled data at the six phases is 
performed, and the Bayesian method is used to correct 
the reconstructed angular distribution results. Six sets 
of corrected data are fitted to obtain six modulation de- 
grees. 


3. Random weights are assigned to the six modulation de- 
grees, with the total sum of the six weights equaling 1. 
The weighted sum yields the modulation degree for this 
sampling. 


After 10,000 samplings, a distribution of the modulation 
degrees is plotted, and a Gaussian fit is applied. The fitted 
sigma represents the Oynfoig- Fig. 17 illustrates the modulation 
distribution of several energy points sampled using the mod- 
ified Bootstrap method from completely polarized data, and 
the Ounfoia obtained from Gaussian fitting. 

For the point 3, the error introduced by parameterization 
can be propagated to the error of the response matrix by pro- 
viding the error of the parameter 7. This error is ultimately 
propagated to the error of the modulation degree. The error 
introduced by parameterization is denoted as para. We ob- 
tain the error of 7 as follows: by making a slight adjustment 
to the value of 7 corresponding to a specific energy point, de- 
noted as Ay. After the adjustment, we incorporate the param- 
eter 7 + A7 into the Star-XP simulation software to simulate 
1,500,000 events, and reconstruct their angular distributions. 
When the x? /ndf between this adjusted angular distribution 
and the angular distribution obtained from the simulation with 
ī equals 1, the Ay represents our estimated error of 7. We 
apply the modified Bootstrap method to resample the data 
10,000 times using the response matrix obtained from the pa- 
rameter value 7 + An, and then perform Gaussian fitting to 
obtain the total error Oynfold+para Or denoted as Gsys. The con- 
tributions of the various error terms at different energy points 
and the modulations corrected are presented in table 1. 

Furthermore, we compared the variation trends of 7 and 
the ratio of position resolution in the Y and X directions 
of the detector at different energies, as shown in Figure 11, 
the trends are in good agreement. Fig.18(a) demonstrates 
the modulations at several energy points before[17] and af- 
ter Bayesian correction, and compares them with the calibra- 
tion results of the IXPE detector[30]. Comparing the modu- 
lations before and after correction, the modulations after cor- 
rection are higher at each energy point than before correc- 
tion. Comparing the corrected modulations with the results 
from IXPE, when the energy is below 4.51 keV, our detector’s 
modulations are lower than the IXPE calibration results. The 
charge induction chip pixel size of IXPE(60um) is smaller 
than the pixel size of our current Topmetal-I(83um), result- 
ing in shorter reconstructed tracks for low-energy photons 
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and higher reconstruction accuracy requirements for resolu- 
tion. Therefore, the modulation of the reconstructed tracks for 
the low-energy part is lower than the IXPE results. However, 
when the energy is higher than 4.51 keV, the corrected mod- 
ulations are higher than the IXPE results, possibly due to the 
better signal-to-noise ratio of Topmetal-II. For longer tracks, 
the pixel resolution no longer plays a decisive role in recon- 
struction accuracy, and factors such as chip noise, the dif- 
fusion coefficient of secondary ionization electrons, detector 
gain, track length, and others begin to have a greater impact 
on the reconstruction. More importantly, after correction, the 
residual modulations of the detector at several energy points 
have been reduced to levels below 1%. In addition, the resid- 
ual modulation result at 5.9 keV in Fig. 18(b) is obtained using 
the etabar at 5.4 keV. The result of the response matrix cor- 
rection at 5.9 keV is 0.24% + 0.59%. The energy resolution 
at 5.4 keV, corresponding to the detector, is approximately 
19.5%[17], and 5.9 keV coincides with the boundary value 
of the 5.4 keV energy resolution. This result indicates that 
the calibration parameter etabar can be extended to the en- 
ergy resolution range of the detector at several energy points, 
while still maintaining good correction results. 
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Fig. 16. The Comparison of the position resolution ratio in the 
Y(ory) and X(orz) directions with 7 at different energies. Blue cir- 
cles represent 7, and green triangles represent ory \Orz- 


V. SUMMARY AND OUTLOOK 


This paper discusses the systematic effects of GMPD and 
corrects the residual modulation of modulation curves caused 
by various systematic effects. GMCP is a prototype detector 
designed for Polar-2/LPD, and the study of GMPD system- 
atic effects is of great significance for the subsequent design 
and performance optimization of LPD, reducing systematic 
effects, and calibrating detector polarization performance. In 
the second section, we list several main systematic effects 
that lead to residual modulation, including differences in gain 
and layout of chip pixels, signal attenuation in electronics, 
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Fig. 17. The modified Bootstrap method is used to sample datasets at different energy points, which are then iteratively corrected through 
Bayesian inference to obtain the weighted modulation degree distribution. The modulation degree distribution is estimated by fitting it with a 
Gaussian function to quantify the error introduced by the Bayesian iteration, denoted as unfold- 


Table 1. The values of 7, Av, Cunfola, para, and the statistical error Ostat at different energy points, as well as the corrected modulations or 
residual of polarized and unpolarized source. 


Energy n Aj} = Ounfold sys Polarization degree Osta  Gtota +Modulation/Residual 
97.4% 0.0079 0.0080 0.2846 
2.98 keV 0.9711 0.00024 0.0009 0.0010 0.0 0.0072 0.0073 0.0075 
99.8% 0.0033 0.0034 0.4743 
4.51 keV 0.9721 0.00068 0.0008 0.0009 0.0 0.0058 0.0059 0.0093 
99.9% 0.0037 0.0039 0.5680 
5.40 keV 0.9819 0.00052 0.0011 0.0013 0.0 0.0057 0.0058 0.0029 
99.8% 0.0038 0.0044 0.6060 
6.40 keV 0.9667 0.00044 0.0016 0.0022 0.0 0.0057 0.0061 0.0084 
99.8% 0.0026 0.0034 0.6117 
8.05keV 0.9619 0.00097 0.0022 0.0023 0.0 0.0039 0.0045 0.0038 
5 045 oF a 3 0.00 5 f t } } f t 
0.30 Pd ł Data after Bayesian correction 
S ł Data before Bayesian correction -0.04 4 
038 ¢ ł IXPE calibration DU1 
i $ scene i i j j à ee ‘au : : 
(a) (b) 


Fig. 18. (a) The red points represent the modulations at different energy points before correction[17], the blue points represent the results 
after correction, and the green points represent the calibration results of the IXPE Detection Unit 1 (DU1). (b) Residual modulation after 
correction of unpolarized data. 
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track truncation, and charge accumulation effects. For these 
known systematic effects, we corrected them through calibra- 
tion, setting threshold conditions, and time positioning. For 
the remaining residual modulation caused by a part of the sys- 
tematic effects, we obtained the response matrix through pa- 
rameterization combined with Monte Carlo simulation and 
used the Bayesian method to eliminate the contribution of 
residual modulation in the modulation curve. The final results 
show that the residual modulation of the data corrected by our 
algorithm has been reduced to below 1% at various calibra- 
tion energy points. The reconstructed modulation degrees of 
data at different polarization phases show good consistency, 
and the polarization degree and modulation degree also ex- 
hibit a good linear relationship. At the same time, we dis- 
cussed the errors of the correction algorithm proposed in this 
paper and compared the corrected modulation results with the 
IXPE calibration results. The data results of GMPD after cor- 
rection by our algorithm show higher polarization detection 
performance than IXPE above 5 keV. 


The results of this paper indicate that the correction algo- 
rithm proposed by us can be well applied to the correction 
of systematic effects in the LPD detector. Additionally, our 
parameterized correction algorithm can naturally be extended 
to the study and correction of oblique incidence systematic 
effects. The correction algorithm that introduces Stokes pa- 
rameters in IXPE is established under the condition of normal 
incidence. When photons are obliquely incident, the descrip- 
tion of photoelectrons using the Stokes parameter system is 
incomplete [31], making it difficult to extend to the correction 
of oblique incidence systematic errors. The large field-of- 
view design of LPD implies that most of the time we need to 
analyze obliquely incident data results. Based on the method 
proposed in this paper, we will also carry out the reconstruc- 
tion and study of oblique incidence systematic effects in the 
future. 
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Appendix A: Bayesian iterative method 


The Bayesian iterative method is a statistical technique 
used to estimate a probability distribution by iteratively up- 
dating prior beliefs with new evidence. It is employed to cor- 
rect for detector effects and estimate the true distribution of a 
physical variable from the measured data. 

In Bayes method, the unfolded distribution, (C4), is given 
by applying the unfolding matrix M;; to the measured distri- 
bution, (£;), as shown in the following equation A1: 


NE 
A(Ci) = X Miyn(E;). (A1) 
j=1 
The unfolding matrix is given by: 
P(E;|C: i 


7 RE P(E;|Ce)no (Ce) 


P(E;|Cx) is the element of response matrix R. € is defined as 
ei = X= P(E;|C;). In the first round of Bayesian iteration, 
no(C) is set based on prior knowledge, while in subsequent 
iterations, no(C) will be replaced by the (C) obtained from 
equation Al. (C) will be updated with each iteration. 

The computation of the error propagation matrix in the 
Bayesian method is given by the following equations A3 A4: 


= Mij + 5 Mixn(Ex) 
k=1 

. 1 Ano(Ci) 
no(Ci) On(E;) 


dn(Ci) 
On(E;) 
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€] ano (Cı) 
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Here, V(n(E;),n(E;)) is computed from the measurement 
data. 
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